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1 Abstract 

In this paper, I propose a realistic interpretation (RI) of quantum mechan- 
ics, that is, an interpretation according to which a particle follows a definite 
path in spacetime. The path is not deterministic but it is rather a ran- 
dom walk. However, the probability of each step of the walk is found to 
depend from some average properties of the particle that can be interpreted 
as its propensity to have a certain macroscopic momentum and energy. The 
proposed interpretation requires spacetime to be discrete. Prediction of 
standard quantum mechanics coincide with predictions of large ensembles 
of particles in the RI. 

2 Introduction 

Despite its general acceptance as a tool to predict the outcome of experi- 
ments with particles and other microscopic objects, standard quantum me- 
chanics (QM) is believed by many to be incomplete or, at least, not fully 
understood. In particular, the "strange" or non-classical phenomena of QM, 
like self-interference and Born's rule, are described in terms of abstract math- 
ematical objects. Although some proponents of the standard or Copenhagen 
intepretation of QM might have been ready to interpret the complex-valued 
wavefunction as a real object [1], wavefunctions are essentially mathematical 
tools serving to calculate probabilities from their square moduli. Contrasting 
to real-valued mathematics and one-to-one mapping between real variables 
and observables of classical theories, the standard description is thus often 
considered as purely operational and thus not realistic. 

Consider the quantum phenomenon of particle self-interference, as illus- 
trated by the double-slit experiment. In the standard picture, self- interference 
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is due to the intrinsic superimposition of the two complex wavefunctions as- 
sociated to the two shts. The wavefunction is beheved to describe every 
single particle, so that a particle is said to pass simultaneously from both 
slits. That implies particle-wave duality, or other hardly imaginable con- 
cepts like diffusion in imaginary time [2]. 

Proponents of alternative views have tried to re-conciliate QM with more 
realistic assumptions concerning particle behavior. A well-known example 
is constituted by hidden-variable theories like De Broglie's pilot wave theory 
or Bohmian mechanics [3], a more recent variant of which is the determin- 
istic trajectory representation [Ij. These theories calculate actual particle 
trajectories [5] that accumulate or rarefact, leading to maxima and minima 
of fringes at the screen behind the slits. However, such an approach lacks 
a plausible mechanism for justifying non-localities, i.e., describing how the 
wave arises and does its guiding. 

The advocates of ensemble or statistical interpretation have recognized 
the standard description of QM as a description that does not apply to 
individual particles but rather to ensembles of similarly prepared particles 
[U [7]. However, ensemble interpretation does not provide any alternative 
way to explain the behavior of single particles, nor it fills the gap between 
the classical domain and the quantum domain 



To explain self-interference in a realistic manner several alternative ap- 
proaches have been developed. A clearly non-exhaustive list includes theo- 
ries involving the influence of the context [HI [10] , nonstandard probability 
definitions [H], nonstandard types of particle trajectories [12], quanta as 
real particles acting as sources of real waves [13], nonlocality and discrete 
spacetime p4]. Despite these efforts, QM is still universally believed to lack 
a convincing realistic interpretation. 

In this paper one possible realistic interpretation is presented. Although 
very simple in its mathematical development and formalism, the proposed 
approach seems to be able to predict self-interference and probability fields 
independently from Schrodinger equation or wavefunctions. Unlike the stan- 
dard picture, only real quantities (actually, integers) are employed. In par- 
ticular, this approach assumes a discrete spacetime. Unlike ensemble in- 
terpretation, it describes single particle trajectories, however, as Markov 
chains, i.e., with an intrinsic randomness. Transition probabilities are simple 
functions of momentum propensity. The latter is randomly determined at 
the preparation phase. Consequently, probability distributions of similarly- 
prepared ensembles of particles are obtained in accord to standard QM pre- 
dictions. 

The paper is organized as follows. Section 1 presents the assumptions 
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concerning spacetime, which naturahy lead to a reahstic interpretation of un- 
certainty principle. Section 2 provides a realistic interpretation of Schrodinger 
equation in terms of Markov chain and transition probabilities. Section 3 
presents the results for a self-interference case, with only one additional as- 
sumption concerning the probability density of momentum propensity as de- 
termined at the sources. Section 4 finally presents a particle-by-particle sim- 
ulation of a double-source experiment, whose a posteriori distribution of ar- 
rival frequency practically tends to coincide with predictions of Schrodinger 
equation. 

3 The lattice: uncertainty principle 

The RI assumes that the spacetime is inherently discrete. Limiting for sim- 
plicity the analysis to one dimension x, the RI assumes that only values 
X = S^X, ^ = 0, 1, 2, . . . and t = tT, r = 0, 1, 2, . . . are meaningful. Noninte- 
ger values of space and time are simply impossible in this picture. The two 
fundamental quantities X and T are the size of the lattice that constitutes 
the space and the fundamental temporal resolution, respectively. 

Under this assumption, a particle trajectory consists of a succession of 
points t} in the spacetime. Advance in time is unidirectional and unitary, 
that is, T + 1 follows necessarily r. Advance in space is still unitary but 
bidirectional. If at a time r a particle resides at the location ^ of the spatial 
lattice, at time t + 1 the particle can only reside at locations ^ + 1, or 
^ — 1. In other words, the local velocity in lattice units f3 = A^/Ar can only 
take the values or —1. 

One must distinguish between local velocity and average or macroscopic 
velocity of the particle. The meaning of the latter is intuitive. Its definition 
in lattice terms is 



The maximum velocity that a particle can reach is the speed of light 
c. Light trajectory in the positive direction corresponds to /3(r) = 1, Vr. 
Consequently to ([T]), one constraint to the fundamental lattice quantities is 
necessarily 

f = c (2) 

Another consequence of ([T]) is that to determine the macroscopic velocity 
of a particle, an observer should wait in principle a time tending to 
infinity. Every observation lasting a finite amount N of time steps will give 





r=l 
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an approximation of v. Consider, e.g., = 1. The observed velocity can 
be +1, or —1. Thus the uncertainty of the macroscopic velocity w is 1 in 
absolute value. For N = 2, the possible results for the observed velocity are 
1, 1/2, 0, —1/2, and —1. Thus the uncertainty of the macroscopic velocity 
is 1/2 in absolute value. Extending these considerations, the uncertainty of 
the macroscopic velocity after an observation lasting time steps is 1/A^ 
in lattice units. 

Moreover, an observation lasting N time steps necessarily implies a 
change in the position of the particle. The span of the particle during the 
observation ranges from NX to —NX. Thus the uncertainty of the position 
of the particle at the end of the observation is obviously 2N in lattice units. 

Using the two results above, and labelling Av{N) and Ax{N) the un- 
certainties of velocity and position as a function of observation time N, the 
relationship 

c 2X2 
Av{N) . AxiN) = -.2NX = ^ (3) 

holds. 

The latter equation resembles the Heisenberg uncertainty principle since 
it fixes an inverse proportionality between the uncertainty with which the 
velocity of a particle can be known and the uncertainty with which its posi- 
tion can be known. Multiplying by the particle mass m, and comparing ^ 
to Heisenberg uncertainty principle, one obtains that the two fundamental 
lattice quantities are related to the Planck constant, 

m-r- = - = — (4) 
T 2 47r ^ ^ 

The term Att (the solid angle of a sphere) is for three dimensional spaces. In 
our example case of a 1-d space, this term reduces to Qi = yJi/2) ' '^'^^re F is 
the Gamma function, i.e., = 2. Thus, combining ([2]) with the accordingly 
modified the values for the fundamental lattice quantities are obtained 
as 

X = ±- (5) 
2mc ^ ' 

and 

^ h 

(6) 



2mc^ 

These values clearly correspond to the Planck length and the Planck 
time, respectively. The role of mass is not completely clear at this point. 
Likely, general relativity might serve to integrate it into the picture. 
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4 Propagation: equation of Schrodinger 



Now I describe the propagation rules of a particle on the lattice in the RI. At 
each time r, the particle might jump to one of the nearest neighboring sites 
of the lattice, or stay at rest. The actual local trajectory is not deterministic, 
i.e., it is not a prescribed function of previous parts of trajectory. Rather, the 
local trajectory has the characteristics of a random walk. This point is very 
important and it implies that an intrinsic randomness affects the particle 
motion. However, in contrast to pure random walks, there is a different 
transition probability for each possible jump. Label the three probabilities 
o, b, and c, respectively. Of course, 

a + b + c = l (7) 

Moreover, the RI assumes that the average or macroscopic momentum, 
p = v/c in lattice units (consider m = 1 henceforth), is imprinted to the 
particle. This imprint is to be attributed to the preparation process. It 
is thus possible to talk about momentum propensity. The quantity p can 
be also reinterpreted as the probability of unitary motion in the positive 
direction, so that 

a — c = p (8) 

Another average or macroscopic quantity is the energy of the particle, 
which in lattice units is written as 



TV 



hm l5^|/?(r)| (9) 



Also this energy can be reinterpreted as the propensity to unitary motion 
in either direction. Thus 

a + c = e (10) 
Combining ([7|)- p0]) . one obtains that 

a=^— , 6=l-e, c=^— (11) 

The energy e must be a function of p. A well-known result of special 
relativity states that energy of a particle is the sum of the rest energy and 
the kinetic energy. In lattice units, that proposition can be written as 

e{p) = (12) 



5 



^=-2 I \ 6 6 

T = T=l 1 = 2 1 = 3 



Figure 1: Particle trajectories after emission from a single source. The graph 
illustrates the neighboring points on the lattice and the three elementary 
probabilities a, b, c. 

Consequently, (jllh can be rewritten as 



step T, determines the trajectory of the particle as a random walk. For 
example, observe Fig. [TJ A particle is emitted from a source located at the 
site ^ = of the lattice with an intrinsic value of p and thus of e, determined 
by the preparation (I will return soon on this point). After one time step, 
r = 1, the particle has a probability a to be at the site ^ = 1, a probability 
b to be at the site ^ = 0, a probability c to be at the site ^ = — 1. After 
two time steps, r = 2, the probabilities for each site from = —2 to = 2 
are as follows: Pp{2) = a^, Pp{l) = 2ab, Pp{0) = 2ac + 6^, pp{—l) = 26c, 
Pp{—2) = (? . Notice that, since the functions a{p) and c{p) are symmetric, 
the probability distribution is symmetric with respect to ^ = 0. Notice 
also that 6^ = 4ac and thus every probability pp is expressed as a unique 
combination of the three elementary probabilities, multiplied by a coefficient. 




(13) 



The equations above, in fixing the probability of each jump at each time 



In general, the recursive expression for the probabihty of finding the 
particle at time r at site ^ is given by 



Pp{^, t) = appiC - 1, r - 1) + bppiC, r - 1) + cpp{^ + 1, r - 1) (14) 

Of course, the start of the recursion is pp{0,0) = 1 for the case of a single 
source. Deriving a closed formula for the probability Pp{C, n) is tedious but 
straightforward at this point. The derivation conduces to binomials and 
binomial coefficients. Instead of attempting to present here the derivation, I 
will show numerical results in the following. However, a special case is easy 
to calculate in closed form. The probability that a particle is at the event 
horizon, i.e., /9p(r, r), is easily calculated as . Similarly, pp{—T,T) = . 

Now, a delicate passage in the theory is introduced. The RI assumes that 
the momentum propensity p is determined randomly during the preparation 
at the source. The probability of releasing a particle with a momentum p is 
uniform over the possible values of p. Since p can vary between -1 and +1, 
its span is 2 and thus the probability density f{p) = 1/2. 

If the source releases a large number of particles in succession, each one 
with a randomly determined value of p, the probability of finding a particle 
at the location ^, t is clearly given by 

Pi(, r) = f{p)Pp{(, r)dp = i ppiC, T)dp (15) 

Let us calculate the integral of (jlSp for the special case mentioned above, 
that is. 



/^(T,r) = -y ^a-dp=-y ^(^^) dp (16) 



Using some elementary mathematics. 



P(7"> t) 



22r+l 



2r + 1 



1 1 1 ^ ^ 

1 22^+1 2r + 1 2r + 1 ^ ' 



This result is not casual. It can be shown numerically that the same proba- 
bility is valid for every other site location comprised between the two event 
horizons, i.e., 

/'(C,^) = ^^,VeG[-r,r] (18) 



It is easily verified that 



Y.p{i,T) = l (19) 
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Now, let us compare this result with the predictions of standard QM, 
i.e., the particular solution of the Schrodinger equation. The wavefunction 
for a free particle is 



1 



ip{k,0)dk 



(20) 



where the wavenumber k is related to the momentum of the particle and 

1 



^(x,0)e-^'=^dx 



(21) 



For a single perfectly localized source at j; = 0, ^{x,0) = 5{x). Thus ([2T 
reads 

ip{k,0) 



1 



and consequently (j20p is rewritten as 

1 



^'(x,t) = 
that, integrated, yields 



exp 



i{kx 



1 



m 



exp 



v/27r V i^t 

The probability density is easily calculated as 

|M'(x,t)|2 



ra 



27rht 



2m 



mx 



m 
hi 



dk 



(22) 



(23) 



(24) 



(25) 



thus it is inversely proportional to time and it does not depend on x. Nor- 
malizing time to lattice units and using @ allows reducing (j25p to 



prob. density = 



(26) 



The probability P(x, t) = \-^{x,t)\'^ X'^ . The result compares with (|17p . with 
2t replacing 2t + 1. The two functions of r are very similar and, indeed, 
practically coincident for r sufficiently large. That could be interpreted in 
the following way: 

The square modulus of the wavefunction predicted by the Schrodinger 
equation is an approximation of the probability in RI. This ap- 
proximation is as better as the measurement is farther from the 
source. 
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Another interesting result arises from the process (|14p. Let us calculate 
the action of the particle, i.e., the energy accumulated by the particle during 
its walk, (t(^, r). Let us introduce also the variable r) = a{^, t)/p{^, t). 
The random process for a is given by 

r) = a(S(C _ i, r - 1) + 1) + r - 1) + c(S(e + 1, r - 1) + 1) (27) 

with the initial condition E(0, 0) = 0. 

For example, one easily obtains a{l,l) = a and thus = 1. The 

action of every particle reaching the point {1, 1} is obviously 1. Generalizing 
this result, clearly a{T,T) = to? and thus S(r, r) = r. For a point like 
{2,0} the prediction is less trivial. The possible values of action can be 2, 
if the particle follows a back-and- forth path, or 0, if it stays at rest for two 
time steps. Using (j27|) . one obtains for this case (t(2,0) = 4ac and thus 
S(2,0) = 2/3, which is a weighted mean between the two possible values of 
action. 

The results easily calculated for other points are listed in Table [H Ob- 
serving the trend of S for a given r, one discovers a clear quadratic depen- 
dency on ^. Indeed, the function E(^, r) is calculated as 

E(e, r) = ^ = Eo(t) + ^ (28) 

) 2t-1 ^ 2t-1 ^ ' 



where Sq is the action at = 0. Equation (p8|) can be easily verified by 
inspection of a few points, provided that {(,,t} > 0. 

Now, compare this result with the phase of the wavefunction (j24p . The 
latter, usually interpreted as the action of the particle is 

S(x.<) = ^ (28) 

which, in lattice units, is 

S{x,t) = 2TT-^ (30) 
2r 

The latter equation corresponds to the second terms in the right-hand side 
of ([28|) . that is, S — Sq, multiplied by 2tt to obtain a phase angle. The 
correspondence is almost perfect, except for the term 2r — 1 that in the 
RI replaces the term 2r predicted by standard QM. For large values of r, 
however, the two results are practically coincident. In other terms: 

The phase predicted by the complex Schrodinger equation is an 
approximation of the action in RI. This approximation is as bet- 
ter as the measurement is farther from the source. 
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Table 1: ( 
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In addition to probability of arrival and action, the source determines 
the whole state of the particle, that is, other observablcs. In order to cal- 
culate a few of them, I find useful to introduce global probability fluxes, 
which are the probability of each possible path ("history", in standard QM 
language) summed over all possible p^s. The flux J^+'°'~^(^, r) describes the 
probability of jumping to site ^ from site ^ — 1, ^, and C + Ij respectively, 
at time r. 

In terms of global probability fluxes, the probability and the action at a 
certain site are also calculated as 

p(e,T) = J^{C,t) + J\C,t) + J-{C,t) (31) 

^(^'^) = ^U+(e,r)[S(e-l,r-l) + l]+ 

+jO(C,r)[S(e,r-l)] + J-(^,r)[S(e + l,r-l) + l]} ^ ' 

Other observablcs of the random walk that can be calculated with this ap- 
proach arc. for instance, the average momentum (the actual observable mo- 
mentum, not the momentum propensity that is a property of the particle) 
at a certain location of the lattice, which is given by 

-f, J+(^,r)- J-(g,r) 

J+(^,r) + JO(^,r) + J-(?,r) ^^^^ 

and is calculated as 

P=- (34) 
r 
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Similarly, the average energy at a certain location is given by 
and it is calculated as 

e = (36) 

r(2T - 1) ^ ' 

Both the latter equation have a counterpart in the variables predicted by 
the Schrodinger equation, in particular in the Bohm interpretation. The 
local average momentum corresponds to the momentum of Bohm hidden 
variable, 

V5 = ^ (37) 

ht ^ ' 

so that the probability flux is j = ;^V5. The local average energy corre- 
sponds to the energy of Bohm hidden variable, that is, time derivative of 
the action 

dS _mx^ _ 1 2 .on^ 

To see the correspondence it is sufficient to write ([37|) ~ (j38|) in lattice units 
or, alternatively, to put m = 1, h = 1 in both of them. As for the other 
observables analyzed above, the correspondence is better as r increases. 



5 Multiple sources: interference 

After having reproduced the predictions of the Schrodinger equation for 
a free particle, I proceed now to a second puzzling aspect of QM: parti- 
cle self-interference. Double-slit experiment usually serves to visualize this 
phenomenon. However, the core of self-interference is isolated and better il- 
lustrated by a double-source preparation. Instead of having a single source, 
a two-slit barrier, and a screen behind the barrier, I represent the same 
process with two independent and mutually alternative sources of particles, 
separated by a certain distance 2dX, and a screen. The sources are equiva- 
lent to very narrow, i.e., punctiform slits. Being in a one-dimensional space, 
the location of the "screen" is clearly fictitious. Pictorially, the geometry 
of the system can be still imagined in two dimensions. One dimension is ^, 
along which the particle move with a momentum propensity p. The second 
dimension is perpendicular to ^ and is traversed by the particle with mo- 
mentum propensity 1 (certainty of advancing in the positive direction). The 
"screen" is thus located at a distance r from the sources. Figure [2] illustrates 
this equivalence. 
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Figure 2: Pictorial equivalence between a double-slit experiment in two 
dimensions and a double-source test case in one dimension. 



The solution of the Schrodinger equation for this case is based on the 
linear superimposition of the two waveforms relative to the two sources, 



^{x,t) 



V2 



(39) 



where both ^i{x,t) and ^2ix,t) are obtained from ()24p by replacing the 
term x"^ (that was valid for a source at x = 0) with a term (x — dX)"^ and 
(x -|- dX)"^, respectively. Thus 



^{x,t) 



1 



m 



exp 



,m{x — dXy 



2ht 



+ exp 



m{x + dXy^ 
' 2ht 



The probability density is given by 



, ^ , , , 9 2m f St — So 
^f(x,t) 2 = —cos^ — 

' ^ ' ^' /it 12 



(40) 



(41) 



where 5*1 and 5*2 are the two independent action values, that is, the phases 
of the two exponentials in (j40p . Finally, the probability in lattice units is 



1 



cos^ vr 



2t 



(42) 



and thus an interference term arises due to the presence of two possible 
sources. The interference is related to the phase difference between the two 
waveforms. 



12 



The representation of the same process in the RI, i.e., in terms of discrete 
spacetime and random walk, would give just the superimposition of two 
probability densities of the type (jlSp , if the same initial conditions are taken 
as in the single-source case. No interference term would arise in this case. 

The key factor to represent self-interference resides in the choice of the 
correct probability density of the momentum propensity p. For a two source 
process, f{p) ^ 1/2. Some values of momentum are more probably than 
others. This fact may seem strange, but actually it is already contained in 
the standard picture of QM [15]. To verify it, it is sufficient to apply (f2T]l 
with ^{x,0) = {5{x - dX) + 5{x + dX))/y/2. The result is 

ifik, 0) = ie-'^^ + e'^'^) (43) 
v27r ^ ^ 



from which the probability function is 

m = l±^^pm (44) 

Recalling that k = p/h, (j44p is transformed in lattice units as 

m = '+°°f""'' (45) 

which is a function oscillating between 1 and with a mean value of 1/2. 

Again, the fact that some momenta are more probable than others should 
not surprise. It is exactly what happens when the two sources are approxi- 
mated by a barrier with two slits at a sufficiently large distance from a single 
source. The two slits filter the momenta spectrum, e.g., favoring the values 
corresponding to the directions of the axes that rely the source to the slits. 
The same behavior is apparently contained in the assumption (j45p . 

Now, simply introducing ()45p into (jlSp . while still using the same process 
(jl4p leading to pp, describes the insurgence of self- interference. This is shown 
numerically in the next section. 



6 Realistic simulation of particle trajectories 

In the last two sections, the predictions of the RI were shown in closed form, 
using mathematical equations in terms of a priori probabilities and proba- 
bility fluxes. The probability of a number of observable were calculated and 
found to be in accord to the predictions of standard quantum mechanics. 
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Table 2: Pseudocode used for the simulations of Fig. \5\ 



i = 

for particle 1 to Np 
i = i + l 

p[i) = random value beween -1 and +1 with prob. den. given by (j45p 
^(0) = random value +1 or -1 with probability 1/2 
s(0) = 
n = 

for time 1 to Nt 
n = n + 1 

V = random value +1, or -1 with prob. given by (jlip 
^(n) =^{n-l)+v 
s{n) = s{n — 1) + \v\ 
end for 

uiaNr)) = y{i{NT)) + 1 
end for 

v{i{NT)) = u{i{NT))/Np 



Now, I will present numerical simulations of the random walk of single par- 
ticles and I will calculate the a posteriori probabilities as frequencies over 
a large number of emissions. Thus, this section is aimed at reproducing 
numerically a true experiment. 

Table [2] shows the pseudocode used for such simulations. The two for- 
cycles are for the successively released Np particles, and for time up to Nt, 
which corresponds to the distance between the sources and the screen in the 
fictitious dimension perpendicular to ^ (see discussion above). Each particle 
experiences the choice of three randomly-selected values: (i) the momentum 
propensity according to /(p), (ii) the source from which it is emitted, and 
(iii) at each time step, its local velocity according to p. The final code line 
represents the counting of the particle that arrive at a certain location at 
time Nt- Prom this number of arrivals, an a posteriori frequency v{^) is 
calculated as the ratio to the total number of particles emitted. 

Pirst, I present the case with only one source at ^ = 0. Pigure [3] shows 
the frequency i^(^) after a time Nt = 300 for different values of Np. As 
the the number of particles emitted in the ensemble increases, a frequency 
distribution builds up. Por large Np, the frequency clearly tends to the a 
priori probability p{^,T = Nt), that is, a constant value given by ([TH]) . 

The case with two sources at ^ = ±1 exhibit a completely different 
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Figure 3: Probability of arrival of a particle emitted at ^ = as a function 
of i after Nt = 300. From top-left to bottom-right, Np = 500, 5000, 10000, 
and 50000, respectively. 



behavior, illustrated by Fig. SI As the the number of particles emitted 
in the ensemble increases, an interference pattern clearly builds up. The 
simulations thus reproduce the outcome of famous experiments like those of 
Merli et al. p,6j and Tonomura et al. [17j . Notice that, despite the fact that 
none of the probability fluxes of the RI can be negative or zero, nevertheless 
there are valleys in the function z^(^), that is, points where the particles 
seldom arrive. 

Figure [5] shows how the frequency v{C) varies with the observation time 
Nt- The number of particles is Np = 50000 in these simulations. The 
curves used as a basis of comparison are the predictions of standard QM, 
namely, the probability function P{^,t = Nt) given by ()42p . Clearly, the 
latter approximate z^(^) as better as the observation time Nt is longer. In 
the case Nt = 300 the function iy(^) is very well smoothed by the analytical 
function P{S,,t = Nt). 
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Figure 4: Probability of arrival of a particle emitted at ^ = ±1 as a function 
of ^ after Nt = 300. From top-left to bottom-right, Np = 500, 5000, 10000, 
and 50000, respectively. 



7 Conclusions 

In contrast to the standard picture of QM, the RI docs not only predict 
the probability of an event over an ensemble of similiarly prepared particles. 
Instead, it also represents a possible behavior of every single particle in the 
ensemble. Probability distribution of the various observables in the state 
are thus obtained a posteriori by evaluating the relative frequency of occur- 
rence of a certain event. Probability distributions can be also determined 
a priori, by applying simple probability conservation equations. Moreover, 
while standard QM calculates the a priori probabilities as squared mod- 
uli of complex-valued wavcfunctions, the RI directly yields the real-valued 
(actually, rational-valued) probabilities, without appealing to mathematical 
abstractions like complex numbers. 

Since no other proposals to describe in a realistic way the beavior of 
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Figure 5: Probability of arrival of a particle emitted at ^ = as a function 
of ^ after four values of Nx- From top-left to bottom-right, Nt = 300, 
100, 40, and 10, respectively. Solid lines: predictions with standard QM, 
P(^,r = Nt). Squares: a posteriori probability z^(^) with the RI. 



single particles exist so far, at least to the knowledge of the author, we 
might be tempted to believe that particles really behave like described in 
the RI. 

The results of this paper have concerned free particles. However, the RI 
seems naturally capable to integrate also external forces into the picture. 
Each interaction of the particle with its sourrounding is indeed expected to 
modify its intrinsic properties, namely, its momentum propensity. 
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